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Abstract 

Using the CLEO II detector operating at the CESR e~^e~ colhder, we have 
measured the structure functions in the decay ir^n^ir^v-r, based on a 

sample corresponding to 4 x 10^ produced r-pair events. We determine the 
integrated structure functions, which depend only on the three pion invari- 
ant mass, as well as the structure functions differential in the Dalitz plot. 
We extract model independent limits on non-axial-vector contributions from 
the measured structure functions as less than 16.6% of the total branching 
fraction, at the 95% confidence level. Separating the non-axial-vector contri- 
butions into scalar and vector contributions, we measure that scalars (vectors) 
contribute with less than 9.4% (7.3%) to the total branching ratio, at the 95% 
confidence level. 

PACS numbers: 13.25.Jx, 13.35.Dx, 14.40.Cs, 14.60.Fg 
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I. INTRODUCTION 



The hadronic structure in the decays tt^tt^tt^i't and tt^tt^tt^i^t- has been the 

focus of several studies in recent years [Q-^. In most of these studies, models are employed to 
characterize the hadronic structure. As data samples have grown, revealing new and complex 
features, the models used have also become more complicated. For example, simple models assume 
the Svr mass spectrum in these decays can be described by a single resonance, the ai(1260) meson. 
DELPHI Q found it necessary to include a radially excited oi meson, the a[ meson. Subsequently, 
CLEO [P, obtained a slightly improved description when including an a'l meson, but found that 
the most significant improvement came about from accounting for a threshold effect in the mass- 
dependence of the oi width due to the opening of the ai K*K decay channel. Similarly, models 
attributing structure in the Svr Dalitz plot distributions solely to decay amplitudes associated with 
ai ptT were found to be insufficient by ARGUS |^ and CLEO Q , who see indications of isospin 
zero contributions, from /2(1270), /o(1370) and /o(400 — 1200) (e.g., a) meson production. The 
models used are not unique, as significant variations in their form and content can lead to similar 
features in the distributions of observable quantities. In addition, no model has so far given a fully 
satisfactory description of the data. 

Thus, although the studies done so far have improved our knowledge, we still do not have a 
clear picture of the hadronic dynamics in r decay to neutrino plus three pions. Furthermore, it is 
difficult to compare results from different experiments and draw additional conclusions since each 
analysis can only be interpreted in the context of the specific model(s) considered. 

J. H. Kiihn and E. Mirkes |^ have constructed a convenient approach for model independent 
studies through the determination of structure functions. In this approach, the kinematics of the 
three pions are studied in a particular reference frame. In this frame, contributions from the 
production of Svr systems with different spin-parity quantum numbers are separated into different 
structure functions, expressed in a completely model independent fashion. This approach also has 
the advantage that the features in the data are condensed into a compact form which can be used to 
test models, alone or in conjunction with data from other experiments. So far, only OPAL has 
determined the structure functions, using the decay — > tt^tt^tt^Vt- They obtained the first fully 
model independent measurement of the signed r neutrino helicity as well as a model independent 
limit on non- axial- vector contributions. 

Here, we present a measurement of the structure functions in the decay — > vr^vr'^vr'^z^T-, based 
on data collected with the CLEO II detector corresponding to 4 x 10^ produced r-pair events. The 
sample of selected events is essentially the same as the one used in our model dependent analysis. 
Consequently, we refer the reader to the article reporting on that work Q for details common to 
both analyses. 

The outline for the rest of the article is as follows. In the next section, we summarize the 
structure function formalism of Kiihn and Mirkes as applied to this decay mode. In Section HI, 
we describe the CLEO II detector and the event selection criteria. The method used to determine 
the structure functions is discussed in Section [V , and the results are presented in Section The 
derived limits on non-axial- vector contributions are given in Section VI. In Section VII we describe 
the sources of systematic error. A summary and discussion of the results follows in Section VIII. 
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II. THE STRUCTURE FUNCTION FORMALISM 



Following the conventions of ref. |^], we denote the momenta of the three pions by q'^ , with 
the momentum of the charged pion given by and the momenta of the two neutral pions by 2 
(1^2 1 > l^iD- The momentum of the total hadronic system is denoted by = Qi + Q2 ^ Q3 its 
invariant mass by \fQ^. 

For the description of the hadronic physics within the context of the structure functions, it is 
convenient, as outlined in ref. Q, to introduce two reference frames y, z) and S'{x' , y' , z'), both 
defined in the three pion rest frame. The two frames are transformed into each other via a rotation 
matrix parametrized by the three Euler angles a, /3, and 7, as illustrated in Fig. |l|. The orientation of 
the reference frame S{x, y, z) is such that the z-axis is parallel to the normal n± = (qi x 52) /\qi x ^2! 
to the three pion decay plane and the x-axis is parallel to ^3, the flight direction of the third pion. 
The reference frame S'{x' ,y' , z') has its z'-axis aligned with hl = —Q/\Q\, i-S- the flight flight 
direction of the laboratory frame as seen from the three pion rest frame. Its azimuthal orientation 
is chosen such that n^, the unobservable flight direction of the r lepton in the three pion rest frame, 
lies in the {x',z') plane. 




In the S{x, y, z) reference frame, using the structure functions Wx and the corresponding lepton 
tensor combinations Lx of Kiihn and Mirkes the differential decay rate for ir^Tr^ir^u can 

be written as follows |^] 

dr{T SiTUr) = C{Q^) ^ Lx{oi, P, 7, 0T, Q'^)Wx{Q^ , si, S2)dQ^ dsids2dad'yd cos /3d cos Or 

X 

with c(q2) = — — cos^ge ^""' 
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and X e {A,B,...,I, SA, SB,..., SG} , 



(1) 



where Gp is the Fermi constant, 6c the Cabibbo angle, and Si the Dahtz plot variables with 
Si = (Qj + Qk)'^ ihJik = l,2,3;i 7^ j 7^ k). The Euler angles a, (3, and 7 as well as the invariant 
mass \/Q^ of the three pion system have been introduced above. The angle 9^ which is the 
angle between the flight direction of the r lepton in the laboratory system and the direction of 
the three pions in the r lepton rest frame, is determined by the energy of the r lepton in the 
laboratory frame, i.e. approximately half the energy of the incoming e~^e~ system, the energy of 
the three pion system in the laboratory frame, and the invariant mass of the three pion system. 
The sixteen structure functions Wx and the corresponding leptonic functions Lx are symmetric 
and antisymmetric combinations of the hadron tensor H^'^ and, respectively, the lepton tensor L^'^ . 
With the exception of the leptonic functions La, Lb, and Lsa all leptonic functions Lx vanish 
after integration over the angles a, /?, 7, and Or- Accordingly, the decay rate in a given bin of 
or {Q'^, 81,82) only depends on the structure functions Wa, Wb, and Wsa- For the details on the 
definitions of the structure functions Wx as well as the leptonic functions Lx the reader is referred 
to ref. II . 

With the exception of the angle a, which determines the azimuthal orientation of the r lepton, 
all variables in Eq. || are measurable quantities in our experiment. Accordingly, in this analysis we 
have to use the differential width dT integrated over a. Performing this integral yields zero for the 
leptonic functions Lsc, Lse, and Lsg (see ref. p). Thus, in our experiment the corresponding 
structure functions Wsc, Wse, and Wsg are not directly measurable. 

Eq. H shows clearly the advantage of the chosen reference frame S{x,y,z). The structure 
functions Wx describing the hadronic physics depend on Q"^ , si, and 82, only. The remaining 
dependence on the angular observables has been rotated into Lx- Accordingly, measurement of the 
set of structure functions Wx{Q^, si, 82) in the reference frame S{x, y, z) yields a model idependent 
determination of the hadronic physics. 

As mentioned above, the structure functions Wx are symmetric and antisymmetric combina- 
tions of the hadron tensor H^^'^ , which on the other hand is derived from the hadronic current /i^ 
by H'^'^ = h^^h*^ . In general /i^ is given by 

= Vl'Fi + V^^F2 + iKf F3 + , (2) 

where the form factors Fi describe the unknown hadronic physics and the vectors V^^ the known 
spin structure. A possible representation is 

XQ{qi-Q3)) 



Vr = q'i-q^-Q^ 



{Q{q2 - qs)) 



Vr = q^ + q^+qi;^Q^. (3) 

The vectors correspond to an axial- vector intermediate hadronic state {e.g. the ai meson), 
to a vector intermediate hadronic state {e.g. the p' meson), and V^^ to a pseudoscalar or scalar 
intermediate hadronic state {e.g. the vr' meson). Accordingly, the form factors Fi^2 describe axial- 
vector contributions, F^ vector contributions, and -F4 pseudoscalar or scalar contributions. As can 
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be seen from Eq. ^, in the reference frame S{x,y, z) the hadronic current /i^ is decomposed into the 
time-Hke component induced by pseudoscalar or scalar contributions, the spatial components 
and /i^ by axial-vector contributions, and the spatial component /i^ by vector contributions. Table 
I summarizes the connection between the different J^-states, i.e. the different components of the 
hadronic current /i^ in the reference frame S{x,y,z), and the structure functions Wx- 



TABLE I. The different J^-states and their connection with the sixteen structure functions Wx- 





J = 


h\ 


JP = 1- 


J = 


WsA 






JP = 1+ 


WsB, Wsc 
WsD, WsE 


Wa 

Wc, Wd, We 




JP = 1- 


WsF, WsG 


Wf, Wg 


Wb 



The vector as well as the pseudoscalar contributions are suppressed due to G-Parity and PCAC. 
Nonetheless they can contribute. For example a possible reaction for the vector contribution in the 
charged pion mode is the decay chain p'^v tt^lou with the uj meson decaying electromag- 

netically to two pions uj — > vr+vr". Ref. estimates that this decay contributes with 0.4% to the 
total TT^TT^TT^u decay rate, in agreement with the measured branching fraction of 0.6% |8|. 

In the neutral mode, — > tt^tt^tt^i', vector contributions can occur via rjir^ mixing: — > -K^ir^rju 
with the r] meson transforming to a pion 7] ^ 2^ ^ n^. The average of the measured branching 
fraction -k^tt^tju from CLEO ^ and ALEPH @ is (0.17±0.03)%. Assuming the ??7r° mixing 

to be of the order of 10~^, ref. Q estimates the vector contribution in the neutral mode to be of 
the order of 10~^. 

States with J = might occur via the vr' intermediate state: ti'^Vt-. Asumming that 

the vr' meson decays subsequently to pvr or frvr, CLEO Q obtained the following upper limits: 

B{t — > Tt'Vt — > pTTVr 'i'KVr) < 1-0 X 10~^ and B{t — > TT'Ur — > (TTTUr — > SlTUr) < 1.9 X 10"^ at the 

90% confidence level. 

Evidently, the sixteen structure functions Wx are not independent from each other. They are 
constructed from the hadronic current h^, which is determined by seven real numbers (an overall 
phase is redundant). Thus, there must exist nine relations among the sixteen structure functions. 
In general these relations are derived from 

H^p ■ H^s = h''h*f^ ■ h^h*^ = h^h*^ ■ h^h*'^ = H^s ■ H^p for all {a, /3, 7, 5} . (4) 

Employing these relations enables one to deduce the structure functions Wsc-, Wse-, and Wsc-, 
despite the fact, as mentioned above, that the leptonic functions Lsc-, Lse, and Lsg vanish after 
integration over the unobservable r azimuth angle a. 



7 



III. DATA SAMPLE AND EVENT SELECTION 



The analysis described here is based on 4 fb^^ of e^e^ colhsion data collected at center-of-mass 
energies 2£^beam of ~ 10.6 GeV, corresponding to 4x 10^ reactions of the type e^e~ — > t^t~ . These 
data were recorded at the Cornell Electron Storage Ring (CESR) with the CLEO II detector |11| 
between 1990 and 1995. 



A. The CLEO II Detector 

CLEO II is a general-purpose large solid angle magnetic spectrometer and calorimeter. Charged 
particle tracking is accomplished with three concentric cylindrical devices: a six-layer straw tube 
array surrounding a beam pipe of radius 3.2 cm that encloses the e'^e~ interaction point (IP), 
followed by two co-axial drift chambers of 10 and 51 sense wire layers respectively. Barrel (| cos 9\ < 
0.81, where 9 is the polar angle relative to the beam axis) and end cap scintillation counters used 
for triggering and time-of-flight measurements surround the tracking chambers. The calorimeter 
comprises 7800 CsI(Tl) crystals, arrayed in projective (toward the IP) and axial geometries in 
barrel and end cap sections respectively. The barrel crystals present 16 radiation lengths to photons 
originating from the IP. 

Identification of iT^TT^ir^i'r decays relies heavily on the segmentation and energy resolution 

of the calorimeter for reconstruction of the vr^'s. The central portion of the barrel calorimeter 
(I cos^l < 0.71) achieves energy and angular resolutions of aE/E{%) = 0.35/^°-^^ -hl.9-0.1^ and 
(70 (mrad) = 2.8/ ^/E + 1.9, with E in GeV, for electromagnetic showers. The angular resolution 
ensures that the two clusters of energy deposited by the photons from a vr'^ decay are resolved over 
the range of vr^ energies typical of the r decay mode studied here. 

The detector elements described above are immersed in a 1.5 Tesla magnetic field provided by 
a superconducting solenoid surrounding the calorimeter. Muon identification is accomplished with 
proportional tubes embedded in the fiux return steel at depths representing 3, 5 and 7 interaction 
lengths of total material penetration at normal incidence. 



B. Identification of Candidate tt^tt tt Ur Decays 

The event selection is nearly identical to that used in our model dependent analysis of the 
-K^-K^T^^v-r decay. Here, we summarize the main features of the event selection procedure. 
For additional details, the article describing the model dependent analysis Q should be consulted. 

To identify events as rr candidates we require the decay of the (denoted as the 'tagging' 
decay) that is recoiling against our signal decay to be classified as Vre^Ve, or Urj-L^fij.- A track is 
identified as an electron if its calorimeter energy to track momentum ratio satisfies 0.85 < E/p < 1.1 
and if its specific ionization in the main drift chamber is not less than 2a below the value expected 
for electrons. It is classified as a muon if the track has penetrated to at least the innermost layer of 
muon chambers at 3 interaction lengths. Thus, we select events that contain two oppositely charged 
barrel tracks (| cos 6*1 < 0.81) with momenta between O.OSE'beam and 0.90£'beam and separated in 
angle by at least 90°, of which one track must be identified as an electron or muon. 

Clusters of energy deposition in the central region of the calorimeter (| cos9\ < 0.71) that are 
not matched with a charged track projection are paired to form ir^ candidates. These showers 
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TABLE II. The four dominant background sources. The total background contribution is 12.6% 
leaving an additional 0.2% background contribution due to multiple sources not listed here. 



fake 








(8.3 ± 0.2)% 


(3.2 ±0.2)% 


(0.5 ±0.1)% 


(0.4 ±0.1)% 



must have energies greater than 50 MeV, and the invariant mass of the photon-pair must lie within 
7.5(T of the TT^ mass where a varies between ~ 4 — 7 MeV. Those vr'^ candidates with energy above 
0.06£'beam after application of a vr'^ mass constraint are associated with any track within 90°. 

A vr^TT^vr^ candidate is formed from a track which has two associated vr'' candidates as defined 
above. If more than one combination of vr" candidates can be assigned to a given track, only one 
combination is chosen: namely, that for which the largest unused barrel photon-like cluster in the 
TT^vr'^Tr'' hemisphere has the least energy. A cluster is defined to be photon-like if it satisfies a 1% 
confidence level cut on the transverse shower profile and lies at least 30 cm away from the nearest 
track projection. 

To ensure that these classifications are consistent with expectations from r decay, events are 
vetoed if any unused photon-like cluster with |cos^| < 0.95 has energy greater than 200 MeV, 
or if any unmatched non-photon-like cluster has energy above 500 MeV. The missing momentum 
as determined using the vr^vr'^Tr'' and tagging systems must point into a high-acceptance region of 
the detector (I cos^jnissl < 0.9), and must have a component transverse to the beam of at least 

0.06£'beam- 

Finally, we define the ttOttO si gnal region to be that where the normalized invariant masses of 
the two photon-pairs, S^^ = (M^^ — ■m^o)/a^^, satisfies —3.0 < < 2.0 for both tt" candidates. 
To estimate the contributions from fake vr'^'s, we also define side and corner band regions using 
—7.5 < S-y-y < —5.0 and 3.0 < S^y < 5.5. The final sample consists of 15849 events in the vr'^vr^ 
signal region. The side and corner band regions contain 1667 and 296 events, respectively. 

The dominant backgrounds are due to mis-identification of r decays to other final states. These 
background modes include (1) decays such as tt^tt^Uj- in which a spurious vr'' is reconstructed 

primarily from secondary clusters arising from interaction of the charged pion in the calorimeter, 
(2) decay modes with three or more vr'^'s, in which the photons associated with the extra vr'^ are 
not detected, (3) the Cabibbo-supressed decay K^ir and (4) the mode — > vr^-fC^z^T-, 

in which the vr^'s originate via the Ks n^n^ decay. The level of contamination from these 
background sources is given in Table 

IV. METHOD 

In our model dependent analysis |Q] , where we were only interested in the relative contributions 
of the amplitudes considered, we used a single entry maximum likelihood fit. Here, we extend the 
single entry maximum likelihood by the normalization to measure the absolute rate in a given bin 
of as well as {Q^, si, 52). As mentioned in section | these rates are determined by the structure 
functions Wa, Wb, and Wsa- 

The likelihood in a given bin j of or (Q^, si, S2) can be written as follows 
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N, 



evt.j 



iLxWxiQ^Si,S2) 



2lnq= V -21nV p"-'''^,',: 



+ 2iV,, (5) 



where N(,^t.j is the number of events in bin and Nj the normahzation of bin j. The upper 
expression in the braces in Eq. |5| is used in the determination of the structure functions Wx 
differential in the Dahtz plot, and the lower expression for the determination of the structure 
functions wx integrated over the Dalitz plot. The leptonic functions Lx are integrated over the 
unobservable r lepton azimuthal angle a including the r pair production P and the factorized 
initial state radiation 

Lx{f3n,0r,Q^) = J hni{K) ■ P ■ Lx{a,P,-f,9r,Q^)ddda, (6) 

where k is the momentum of the radiated photon. For the initial state photon spectrum we 



use the formula obtained in ref. [12|. Overall constant factors of the r pair production have been 
factored out and, accordingly, P denotes only the functional form of the r pair production with 
P = (1 + cos^ ^T-^)/2, where Orr is the polar angle of the r lepton with respect to the beam axis. 
The integrals over a and k in Eq. ^ are evaluated numerically. 

We have included the initial state bremsstrahlung in the leptonic functions Lx to account for 
the fact that the r energy Er in the laboratory frame, needed for the evaluation of the angle Or 
(Eq. HI), is in general not given by Ej- = but by Ej- = ^/s'/2, where \/s' is the center of mass 

energy after radiation. 

The inclusion of the r pair production corrects for the fact that the (1 + cos^Orr) distribution 
of the r leptons induces an anisotropic distribution of the r momenta that lie on a cone around the 
three pion momentum. Azimuth angles that result in r momenta closer to the beam are preferred 
over azimuth angles further away from the beam. In the absence of any selection cuts this anisotropy 
vanishes after integration over the orientation of the r lepton or, equivalently, the orientation of 
the three pion system in the laboratory frame. However, the selection cuts applied yield a residual 
anisotropy affecting the leptonic functions Lx- 

Monte Carlo studies showed that the neglect of the r pair production as well as of the initial 
state radiation results in a bias in the determination of the structure functions. In both cases the 
bias obtained is around 2%. Although the statistical uncertainty in our results is much larger than 
the effect of r pair production, we expect that future analyses will not be limited by statistics and 
that pair production dynamics will be more important. 

For a given bin j in or in (Q^, si, S2) the normalization Nj in Eq. ^ is chosen to be 

iV, = -^xr, = -^xC(Q)Ey| ^^^^^Q2^ jd,dcosPdcoser, (7) 

where 7, Or, and C{Q'^) have been introduced above (see Eq. ||). The total number of events 
is denoted by N^vt- The partial width, F, of the lepton decaying into vr^Tr^vr'^^'T- is given by 
the total width Ttot = ^/^r times the branching fraction B{t^ — > Tr^T^^ir^Vr)- We use the world 
average values of Fto* = (2.27 ± 0.01) x 10^3 g ^(^t ^ ^TyrOvrOi/^) = (9.15 ± 0.15)% 
to determine the partial width F. The factor /^e/ = 0.996 corrects for events that are outside of 
our selected region of 0.5 GeVVc^ < < 2.75 GeVVc^ i.e. Y.j = fsei x F. 

As in our model dependent anaylysis Q, we account for the four main background sources, 
listed in Table by extending the likelihood of Eq. |5| as follows 
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C (1 ^f^O ^iiT ^K-KTT (^KsTT^ ^signal 

where Csignai is the hkehhood of the signal events, J~-f ^ the likelihood of the fake background, 

the likelihood of the 7r^(37r'^)^'T- background, Ck-k-jt the likelihood of the K 

background, and Ckstt the likelihood of the r"^ Kst^^Vt background. The corresponding back- 
ground fractions aj are given by Table ^ The likelihood >C/ ^ of the fake pion background is 
approximated by the Dalitz plot distributions of events populating the t:^ mass side bands. We 
model the likelihood C^t^ of the four pion background with the decay /9(1450) — > pa using a S'-wave 
amplitude. To estimate the systematic error arsising from the large theoretical and experimental 
uncertainites of the four pion matrix element we also consider the decay /?(1450) — > aivr (5-wave). 
The integration over the lost vr'^ meson, needed for the evaluation of C^-,^, is done numerically 
taking into account efficiency. The background K^n^TT^Ur is modeled by the decay chain 

— > K^u-r, K*^TT^ (S'-wave), where the Ki meson is parametrized by a superposition 

of the ifi(1270) and i^i(1400) Breit Wigner functions. Finally, the Kstt^i't background is 

parametrized by the decay chain K*^Vt-, K*^ — > Kgir^ (P-wave). The mass distribution 

for the Ks — > tt^tt^ decay is parametrized by a Gaussian, where the mean and the width are taken 
from data. 

To determine the structure function integrated over the Dalitz plot we subdivide our data sample 
in nine equidistant bins of Q^. The considered range is 0.50 GeV^/c^ < < 2.75 GeV^/c^ 
yielding a bin width of 0.25 GeV^/c^. In each of the nine bins we fit for the integrated structure 
functions wa, wc^ '^Di ^-nd we- In this case non- axial- vector contributions are not taken into 
consideration, since, due to the loss of information by integrating over the Dalitz plot, we are 
unable to resolve these contributions. 

Table III shows in detail the binning in (Q^, si, S2) used to extract the structure functions Wx 



differential in the Dalitz plot. In total we have 34 bins. Instead of fitting directly for the structure 
functions Wx and employing explicitly the nine relations among the structure functions (see Eq. ^) , 
we introduce in each bin of {Q'^, 81,82) seven fit parameters: the real 3ft(/io) and imaginary part 
^(/lo) of the time-like component of the hadronic current h^, the real part K(/ii) of the x-component, 
the real 3f?(/i2) and imaginary part 9(/i2) of the y-component, and the real ^(/is) and imaginary 
part ^(/is) of the z-component. The imaginary part of the x-component 9(/ii) is chosen to be 
zero. From the measured hadronic current /i^ in a given bin of {Q'^ , 81, 82) we then calculate the 
structure functions Wx taking into account the full covariance matrix for the error evaluation of 
the structure functions Wx- 



V. MEASUREMENT 

The fit results obtained for the integrated structure functions wa, wc-, wd, and w^; in compari- 
son with the model obtained in our model dependent analysis Q as well as the model of J. H. Kiihn 
and A. Santamaria |jl^ (KS) are shown in Fig. ^. The observed steep falling of the structure func- 
tion WA around = 2 GeV^/c'^, which might be an indication of the K*K threshold opening, is 



well reproduced by the fit obtained in our model dependent analysis Q. The KS-model [13|, which 
does not include the K*K threshold in the total width of the ai meson Breit Wigner function, falls 
less steeply around = 2 GeV^/c'*, as can be seen from Fig. |2|. Owing to the large statistical 
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TABLE III. Employed binn 
The Dalitz plane is symmetrize 
space boundaries (see e.g. ref. | 


ing in {Q'^ , si, S2) as used in the fits differential in the Dalitz plot, 
d with si > S2. Depending on si, 32^°"^ is either given by the phase 
14 , section 35.4.3.1) or by the condition si > S2. 


Bin 


Q2 [GeV^] 


si [GeV^] 


S2 [GeV^] 




1 

2 


0.50-0.75 


0.00 - 0.35 
0.35 - 0.75 


0.00 - s^"^ 
0.00 - s^"^ 


192 
204 


3 
4 
5 
6 


0.75 - 1.00 


0.00 - 0.48 
0.48 - 0.56 
0.56 - 1.00 
0.32 - 0.68 


0.00 - 0.32 
0.00 - 0.32 
0.00 - 0.32 
0.32 - s^'^^ 


527 
575 
492 
512 


7 
8 
9 
10 
11 
12 


1.00 - 1.25 


0.00 - 0.53 
0.53 - 0.59 
0.59 - 0.66 
0.66 - 1.25 
0.40 - 0.56 
0.56 - 0.85 


0.00 - 0.40 
0.00 - 0.40 
0.00 - 0.40 
0.00 - 0.40 
0.40 - s^""^ 
0.40 - s^""^ 


680 
619 
663 
614 
645 
667 


13 
14 
15 
16 
17 
18 


1.25 - 1.50 


0.00 - 0.56 
0.56 - 0.64 
0.64 - 0.74 
0.74-1.50 
0.52 - 0.66 
0.66 - 0.98 


0.00 - 0.52 
0.00-0.52 
0.00-0.52 
0.00-0.52 
0.52 - s^"^ 
0.52 - s^'^^' 


667 
675 
643 
663 
665 
640 


19 
20 
21 
22 
23 
24 


1.50 - 1.75 


0.00 - 0.61 
0.61 - 0.73 
0.73 - 0.89 
0.89 - 1.75 
0.57-0.77 
0.77- 1.18 


0.00-0.57 
0.00 - 0.57 
0.00 - 0.57 
0.00-0.57 
0.57 - s^"^' 
0.57 - 3^^"= 


496 
481 
442 
440 
490 
498 


26 
27 
28 


1.75 - 2.00 


n nn n 79 
0.72 - 1.00 
1.00 - 2.00 
0.66 - 1.34 


n on n fifi 

u.uu — u.uu 
0.00 - 0.66 
0.00 - 0.66 
0.66 - s^"^' 


379 
360 
373 


29 
30 


2.00 - 2.25 


0.00 - 1.05 
1.05 - 2.25 


0.00 - 3^"^ 
0.00 - s^''^ 


320 
300 


31 
32 


2.25 - 2.50 


0.00 - 1.25 
1.25 - 2.50 


0.00 - s^"^ 
0.00 - s^''^ 


109 
124 


33 
34 


2.50 - 2.75 


0.00 - 1.20 
1.20 - 2.75 


0.00 - s^"^ 
0.00 - s^''^ 


55 
47 



errors obtained on the structure functions wc, wd, and we the comparison of the fit results and 
the models is not very conclusive for those three functions. 

The fit results for the axial- vector part of the hadronic current (Table |) in bins of (Q^, 31,32) 
are shown in Fig. The measured non-axial-vector contributions of the hadronic current are 
plotted in Fig. |^. 

As mentioned in section IV, we derive the structure function Wx differential in the Dalitz plot 
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FIG. 2. The integrated structure functions wa, wc, wd, and we- The filled points represent 
the data. The smaller error bars indicate the statistical errors. The distance between the smaller 



and larger error bars shows the systematic errors (see Section VII ). The solid line is the model as 
obtained in our model dependent analysis, whereas the dotted line corresponds to the KS-model. 



from our measured hadronic current including the full covariance matrix. The results obtained on 
the axial- vector induced structure functions Wa, Wc, Wd, and We are shown in Fig. ^, whereas 
Fig. ^ shows the non- axial- vector structure functions Wb, Wsa, Wse, and Wsg- The structure 
functions that have their origin in the interference between the vector and axial- vector contributions 
of the hadronic current as well as the ones due to the interference between the scalar and axial- vector 
components of the hadronic current are plotted in Figs. ^ and |8|, respectively. 

The KS-model H, as can be seen from Fig. |3|, fails to describe the details of the hadronic 
current components and 3f?(/i2). This discrepancy between the data and the KS-model |13| 

is also reflected in the functions Wa and Wc, shown in Fig. 

As illustrated by Fig. ^ and, correspondingly, by the Figs. ^, 0, and ^, we do not measure 
significant scalar or vector contributions. The null hypotheses for scalar and vector contributions 
as well as for scalar and vector contribution combined yield in all three cases significances below 
three standard deviations. 

Tables summarizing our results on the structure functions as well as on the hadronic current 
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components are available on the WWWQ. 
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FIG. 3. The real and imaginary parts of the axial-vector induced components of the hadronic 
current as measured in our fits differential in the Dalitz plot. For the definition of the bin number 



the reader is referred to Table III. The filled points represent the data. The smaller error bars 
indicate the statistical errors. The distance between the smaller and larger error bars shows the 



systematic errors (see Section VII). The solid line is the model as obtained in our model dependent 
analysis, whereas the dotted line corresponds to the KS-model. 



VI. LIMITS ON NON AXIAL- VECTOR CONTRIBUTION 



The fit results on the scalar components of the hadronic current, i.e. the time-like components 
5R(/io) and 9(/io), yield the following branching fraction for scalar contributions 



5(rT ^ TT^FTrOvrOzy) 



(5.2 ±2.5 ±0.3)% 



(9) 



^See |http:/7 wwvy. Ins. Cornell. edu/pubIic/CLEO/analysis/results/tau-struct/|. 
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FIG. 4. The real and imaginary parts of the non-axial-vector induced components of the 
hadronic current as measured in our fits differential in the Dalitz plot. For further explanations see 
comments in caption for Fig. H. 



where the first error is statistical and the second error systematic. The individual contributions to 
the systematic error are: ±0.29% due to efficiency, ±0.06% due to Monte Carlo statistics, ±0.03% 
due to the energy scale of the photons, and ±0.01% due to background. With this we obtain the 
following upper limit 

^Kr/' Z^'Il^'^y < 9-4% 95% CL. (10) 

Using the measured contributions of the vector components ^(h-^) and ^(/la) we measure the 
branching ratio for vector contributions to be 

^00.^ = 4.2 ± 1.9 ± 0.2 % , 11 
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FIG. 5. The axial-vector induced structure functions Wa, Wc, Wd, and We as obtained from 
our fits to the real and imaginary parts of the hadronic current. For further explanations see 
comments in caption for Fig. H. 



where the first error is statistical and the second error systematic. The systematic error is given 
by the quadratic sum of a ±0.15% contribution due to efficiency, a ±0.10% contribution due to 
Monte Carlo statistics, a ±0.02% contribution due to the energy scale of the photons, and a ±0.01% 
contribution due to background. The upper limit on vector contributions is 

/ ^ n n ^ < 7-3% at 95% CL. (12) 

Combining the results on the scalar and vector components of the hadronic current we obtain 
the following branching ratio for non-axial- vector contributions 

0{t+ — > -K+T^^TT^V) 
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FIG. 6. The non-axial- vector induced structure functions Wb, Wsa, Wsf, and Wsg as obtained 
from our fits to the real and imaginary parts of the hadronic current. For further explanations see 
comments in caption for Fig. H. 



where the first error is statistical including the correlations between the vector and scalar contri- 
butions and the second error systematic. The systematic error is given by the quadratic sum of the 
systematic error of ±0.30% from the scalar contribution and the systematic error of ±0.18% from 
the vector contribution. The upper limit on non-axial-vector contributions is 

B(t^ Non-Axial u tt^tt'^tt^v) ^ . ^ 

^ — ^ n n , < 16.6% at 95% CL. 14 

The corresponding upper limit obtained by the OPAL collaboration Q is 



< 26.1% at 95% CL. (15) 
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FIG. 7. The structure functions Wf, Wg, Wh, and Wi — induced by the interference between 
the vector and axial-vector components of the hadronic current — as obtained from our fits to the 
real and imaginary parts of the hadronic current. For further explanations see comments in caption 
for Fig. 0. 



VII. SYSTEMATIC ERRORS 

We consider the following sources of systematic errors: 

• Overall Normalization: As mentioned in section we normalize our likelihood to the world 
average branching ratio of B{t^ ir^ir^TT^UT-) = (9.15 it 0.15)%. The error on the branching 
ratio results in an overall systematic error of our normalization. In addition, the extrapolation 
of our measurement into the regions 

Q2 < 0.5GeVVc^ and > 2.75GeVVc^ i.e. the 
uncertainty on fsei in Eq. 0, also yields an error on our normalization. We conservatively 
estimated this error to be as large as the error on the branching ratio. Accordingly, we 
have an overall systematic error on all structure functions of 3% and, in case of the hadronic 
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FIG. 8. The structure functions Wsb-, Wsc-, Wsd, and Wse — induced by the interference 
between the scalar and axial-vector components of the hadronic current — as obtained from our 
fits to the real and imaginary parts of the hadronic current. For further explanations see comments 
in caption for Fig. |3[ 



current components, an overall systematic error of 1.7%. 

Monte Carlo Statistics: The limited statistics in the calculation of the normalization integrals 
of the likelihood function for the different bins in as well as in (Q^,si,S2) results in a 
systematic error of our measurement. We estimated this error by subdividing our Monte Carlo 
sample, which has approximately ten times the size of our data sample, in six independent 
sub-samples. The variation in the fit results obtained are used as an estimate for this error. 

Efficiency: The not fully accurate detector simulation of the efficiency of the momenta of the 
three pions and the angle between the pions as well as the invariant mass of the three pions 
can also effect our measurement. We estimate this uncertainty by reweighting our events as 
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a function of the momenta of the three pions, the angle between the pions, and the invariant 
mass of the three pions. Since our hkehhood is normahzed to the world average branching 
ratio Q, the error in the overall efficiency does not affect our measurement. Accordingly, 
we only account for relative changes of the efficiency. 

• Background: The amount of background as well as the background models used in the 
likelihood have uncertainties. The systematic error associated with these uncertainties have 
been estimated by varying in the likelihood function the background contribution and, in 
case of the four pion background, also the model, i.e. (1) /9(1450) pa, (2) yo(1450) — > aivr, 
and (3) combinations thereof (see section 

• Photon Energy Scale: The energy scale of the photons that form the neutral pions has an 
uncertainty of around 0.3%. The resulting systematic uncertainty in our measurement has 
been estimated by rescaling the energies of the neutral pions by 0.3% and refitting our data. 

In addition to the systematic uncertainties listed above, we also have an uncertainty due to the 
finite detector resolution. Refitting our events including the error matrix on the charged track and 
the photons (scaled by a factor of four) showed that this source is negligible compared to the other 
sources and has, therefore, not been taken into account. This coincides with our findings of our 
model dependent analysis 

In general all five sources result in systematic errors that have roughly the same order of 
magnitude. However, the errors due to the overall normalization and the Monte Carlo statistics 
yield the largest systematic uncertainy for most of the bins. 

The systematic errors obtained are shown in Figs. ^, ^, ^, ^, ^, ^, and ^ as extension to the 
error bars. 

VIII. SUMMARY 

We have measured the integrated structure functions wa, wc-, wd, and we- In addition, we 
determine the seven non-trivial real and imaginary components of the hadronic current in bins 
of {Q"^, 81,82). The results obtained on the hadronic current enable us to determine all sixteen 
structure functions Wx differential in the Dalitz plot. No significant non-axial- vector contributions 
have been found. 
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